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Abstract 

To study discrete dynamical systems of different types — deterministic, statistical and 
quantum — we develop various approaches. We introduce the concept of a system of 
discrete relations on an abstract simplicial complex and develop algorithms for anal- 
ysis of compatibility and construction of canonical decompositions of such systems. 
To illustrate these techniques we describe their application to some cellular automata. 
Much attention is paid to study symmetries of the systems. In the case of determin- 
istic systems, we reveal some important relations between symmetries and dynamics. 
We demonstrate that moving soliton-like structures arise inevitably in deterministic 
dynamical system whose symmetry group splits the set of states into a finite number 
of group orbits. We develop algorithms and programs exploiting discrete symmetries 
to study microcanonical ensembles and search phase transitions in mesoscopic lattice 
models. We propose an approach to quantization of discrete systems based on intro- 
duction of gauge connection with values in unitary representations of finite groups — 
the elements of the connection are interpreted as amplitudes of quantum transitions. 
We discuss properties of a quantum description of finite systems. In particular, we 
demonstrate that a finite quantum system can be embedded into a larger classical sys- 
tem. Computer algebra and computational group theory methods were useful tools in 
our study. 
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1. Introduction 

There are many reasons — physical, mathematical, and conceptual — to study discrete 
structures. Discrete systems are important in applications — nanostructures, for exam- 
ple, by their nature are discrete, not continuous, formations. From a fundamental point 
of view, there are many philosophical and physical arguments that discreteness better de- 
scribes physic^] at small distances than continuity which arises only as approximation or 
as a logical limit in considering large collections of discrete structures. As a recent de- 
velopment, let us mention much-discussed E. Verlinde's thermodynamic (entropic) deriva- 
tion HI of gravity and Newton's law of inertia from G. 't Hooft's holographic principle 
||2). The holographic principle conjectures that it is possible to describe physical events in 
a three-dimensional volume fully by a theory on its temporally varying two-dimensional 
boundary — holographic screen — containing finite number of discrete degrees of free- 
dom. Entropy of these degrees of freedom, i.e., number of bits N, is proportional to the 

area A of the screen: N = ^^-Jl In more speculative sense, the whole universe is a finite 
two-dimensional information structure on the cosmological horizon, and observable three 
dimensions are only an effective description at macroscopic scales and at low energies. 
Verlinde shows that the laws of Newton and the Einstein equations come out directly and 
unavoidably from the holographic principle. The gravity appears to be an entropic force 
arising in systems with many degrees of freedom by the statistical tendency to increase its 
entropy — like osmosis or elasticity of polymers. Verlinde derived his results combining 

holography [n = , the equipartition rule (assumption on even distribution of energy 

over N bits), 1st law of thermodynamics (dE = TdS — Fdx) and several additional stan- 
dard relations. To introduce thermodynamics, i.e., to construct canonical partition function, 
there is no need to know details of microscopic dynamics. It suffices to know about energy 
and number of states. Of course, the fundamental problem about laws governing bit dy- 
namics on holographic screens remains unsolved. Since Planck scales are experimentally 
unavailable — the Planck length is about 10 -35 meters, i.e., far below the spacial resolu- 
tion of particle accelerators (nowadays about 10~ 18 meters) — the construction and study 
of various discrete dynamical models is one of the possible approaches. 

In this chapter we consider three types of discrete dynamical systems: deterministic, 
mesoscopic statistical and quantum. 

We begin with a general discussion of dyscrete dynamical systems. The most funda- 
mental concepts are a discrete time and a set of states evolving in the time. A space is 
considered as a derived concept providing the set of states with the specific structure of a 
set of functions on the points of space with values in some set of local states. We give an 
illustration of how a space-time may arise in simple models of discrete dynamics. Then we 
discuss symmetries of space and local states and how these symmetries can be combined 
into a single group of symmetries of the system as a whole. 

We introduce the concept of a system of discrete relations on an abstract simplicial 

'Of course, the question of "whether the real world is discrete or continuous" and even "finite or infinite" 
is rather metaphysical, i.e., neither empirical observations nor logical arguments can validate one of the two 
adoptions — this is a matter of belief or taste. 

2 In theories with emergent space this relation may be used as definition of area: each fundamental bit 
occupies by definition one unit of area. 
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complex EJIH, and explain how any system of discrete relations — subsets of Cartesian 
products of finite sets — acquires the structure of an abstract simplicial complex. This 
general concept covers many discrete mathematical structures. In particular, it can be con- 
sidered as generalization of cellular automata or as a set-theoretical analog of systems of 
polynomial equations — if all factors of the Cartesian product are sets with the same num- 
ber of elements and this number is prime power, than any relation can be expressed by 
polynomial equation. We describe algorithms for analysing compatibility and constructing 
canonical decompositions of discrete relations. As an illustration, we give results of ap- 
plication of the algorithms to some cellular automata, namely, Conway's automaton Game 
of Life and Wolfram's elementary cellular automata. For many of the latter automata the 
canonical decomposition allows to obtain either general solutions in closed form or impor- 
tant information on their global behavior. 

Symmetry is a property of fundamental importance for any mathematical or physical 
structure. Many real world discrete systems, e.g., carbon nanostructures like graphenes and 
fullerenes, are highly symmetric formations. Symmetries play essential role in the dynam- 
ics of the systems. In this chapter we consider connection between symmetries of discrete 
dynamical systems on graphs — 1 -dimensional simplicial complexes — and their dynamics 
OHl. In the case of deterministic dynamical systems, such as cellular automata, non-trivial 
connections between the lattice symmetries and dynamics are revealed. In particular, we 
show that formation of moving soliton-like structures — typical examples are "spaceships" 
in cellular automata — is a direct result of the existence of non-trivial symmetry. 

We developed also algorithms exploiting symmetries for computing microcanonical 
partition functions and for searching phase transitions in mesoscopic lattice models. 

We consider a class of discrete dynamical models allowing quantum description (71. 
Our approach to quantization consists in introduction of gauge connection with values in 
unitary representation (not necessarily 1 -dimensional) of some group of internal symmetries 
— the elements of the connection are interpreted as amplitudes of quantum transitions. The 
standard quantization is a special case of this construction — Feynman's path amplitude 
e tfLdt can k e interpreted as parallel transport with values in (1 -dimensional) fundamental 
representation U(l) of the group of phase transformations. For discrete systems it is natural 
to take a finite group as the quantizing group, in this case all manipulations — in contrast 
to the standard quantization — remain within the framework of constructive discrete math- 
ematics requiring no more than the ring of algebraic integers (and sometimes the quotient 
field of this ring). On the other hand, the standard quantization can be approximated by 
taking 1 -dimensional representations of large enough finite groups. 

Any approach to quantization leads ultimately to unitary operators acting on a Hilbert 
space. We discuss peculiarities of quantum description of finite systems, under the assump- 
tion that the operators describing quantum behavior are elements of unitary representations 
of finite groups. We show that in this case any quantum problem can be embedded into a 
classical one with a larger space of representation. 

Computer algebra and computational group theory [8] methods turned out to be quite 
useful tools in our study of discrete systems. 
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2. Discrete Dynamics 

Generally, discrete dynamical system is a set S = {si, . . . , sjy 5 } of distinguishable states 
evolving in discrete time t G T = Z = {. . . , —1, 0, 1, . . .}, i.e., evolution or history is an 
element of the set £ = S T . Dynamics is determined by some evolution rule connecting 
the current state s t € S of the system with its prehistory st-i, s t ~2, st-3, ■ ■ ■ Different 
types of evolution rules are possible. We shall consider here the following types of discrete 
dynamics. 

• Evolution rule of deterministic dynamical system is afunctional relation. This means 
that the current state is a function of the prehistory: 

S t =F(s t -l, S t -2,S t -3, ...)• (!) 
Cellular automaton is a typical example of deterministic dynamical system. 

• Statistical lattice model is a sort of non-deterministic dynamical system. This is a 
special case of Markov chain. In statistical lattice model transition from one state to 
any other is possible with probability controlled by a Hamiltonian. 

• Quantum system is another important type of non-deterministic dynamical system. 
The probabilities of transitions between states are expressed in terms of complex- 
valued transition amplitudes. 

Symmetries play an important — central in the case of quantum systems — role in 
dynamical systems. So we assume the existence of a non-trivial group W = {wi = 1, W2, 
. . . , wat w ]|j acting on the set of states S: W < Sym (S). Action of the group W splits the 
set of states S into orbits of different sizes: S = |J 0\ (disjoint union). 

i 

2.1. Discrete Dynamical Models with Space 

In applications the set of states S usually has a special structure of a set of functions on 
some space.The following constructions form the basis for all types of dynamical systems 
we consider in this chapter: 

1. Space is a discrete (basically finite) set of points X = {xi,X2, . . . ,xjv x } provided 
with the structure of an abstract regular (/c-valent) graph. 

2. Space symmetry group G = {gi = l,g2, • • ■ >g7V G } * s tne g ra ph automorphism 
group: G = Aut(X) < Sym (X). We assume that G acts transitively on X. 

3. Local space symmetry group is defined as the stabilizer of a vertex jq in the space 
group G: <? G Gi oc = Stabc (x) means Xj<7 = Xj@ Due to the transitivity all such 
subgroups are isomorphic and we shall denote the isomoiphism class by Gi oc . This 
is subgroup of the space symmetry group: Gi oc < G. 

3 We denote the identity elements by 1 for all groups throughout this chapter. 

4 We write group actions on the right. This, more intuitive, convention is adopted in both GAP and Magma 
- the most widespread computer algebra systems with advanced facilities for computational group theory. 
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4. Points x G X take values in a finite set £ = {a±, 02, . . . , <tn s } of Zoca/ states. 

5. Internal symmetry group T = {71 = 1, 72, . . . , 77v r } is a group T < Sym (£) acting 
on the set of local states S. 

6. States of the whole system are functions <r(x) £ £ x = 5, and the set of evolutions 
takes the form £ = (S X ) T = S XxT . 

7. We define the whole symmetry groups W unifying space G and internal F symmetries 
as equivalence classes of split group extensions of the form 

1^T X ^W^G^1, 

where T x is the set of T-valued functions on X. (More detailed description of this 
construction see in Sect. 12.1.3.1 ) 

The separation of the set S into "space" and "local states" is not fundamental — it is model- 
and interpretation-dependent. An example of a system with a somewhat non-standard no- 
tion of space is a quantum computer. Here the space X is the set of Nx qubits, the set of 
local states S is {0, 1}. The whole set of states S = {0, 1} X contains 2 JVx elements. 

2.1.1. Example of Discrete Model with Emergent Space-time. 

Modern fundamental theories, in particular the string theory, provide evidence that space 
is an emergent phenomenon Q, arising from more basic concepts. We demonstrate here 
that if we have a concept of time then discrete space-time structures may arise under very 
simple and general assumptions. It is sufficient to have a time-labelled sequence of events 
and ability to distinguish different types of the events. Then space dimensions arise as the 
counters of events of different types. 

Let us consider a set of states (symbols) £ = {a\, 0%, . . . , o~n+i} and assume that it is 
possible to observe the sequences (histories) h = so, s\ ■ ■ ■ St, where Sj 6 E. Let us define 
a space-time point p as equivalence class of sequences with equal numbers of occurrences 
of each symbol, i.e., p is a "commutative monomial" of the total degree t described by 
N + 1 non-negative integers: p = (m, . . . , njy+i) , n\ + • • • + n^v+i = t, n» £ Z>o is 
multiplicity of symbol cjj in the history h. The concepts of "causality" and "light cones" 
arises naturally. The "speed of light limitation" is simply impossibility to get more than t 
symbols ("perceptions") in t observations — in terms of monomials the "past light cone" is 
the set of divisors of the monomial p, the "future light cone" is the set of its multiples, see 

Fig.m 

The union of all possible histories form a causal network. As to modelling continu- 
ous Euclidean spaces by this structure, the system of discrete points can be embedded into 
a continuum in many different ways: as a set of discrete points into a continuous space 
of arbitrary non-zero dimension^, as a network into a three-dimensional space^. To sepa- 
rate space from the space-time one should introduce a rule identifying points at different 

5 For example, the map p — > aini + - ■ - + Qjv+injv+i £ R 1 , where are independent irracionals, provides 
one-to-one embedding of the set of points into R 1 . 

6 A network, as a locally finite 1 -dimensional simplicial complex, can always be embedded into R 3 . 
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times. The identification can be any causality-respecting projection onto the set t = const. 
To construct an illustrative discrete model of this section we use the following projection. 
Let us identify the symbols o~i with N + 1 unit vectors forming regular simplex in an 
iV-dimensional Euclidean space. These systems of vectors (network generating sets) look 



like 



for N = 1,2,3, respectively. The space lattices generated 



by these sets in four time steps (t = 4) for the cases N = 2 and N = 3 are shown in the 
figure 





With these prerequisites, let us construct a simple physical model in 1 -dimensional 
space (N = 1). We have E = {0-1,(72} = {— >, <— }, t = n\ + rii € Z>o- Let us add a 
little physics by imposing the structure of Bernoulli trials on the sequences h = sqSi ■ ■ ■ St. 
Namely, let us introduce probabilities p\ and P2 (j>2 + Pi = 1) for possible outcomes o\ 
and o"2 of a single trial. The probability of a separate history h is described by the binomial 
distribution 

P{n l ,n 2 ) = j — 7—p l 1 p 2 . (2) 

ni\n2\ 

From this model we can see that the behavior of a discrete system may differ essentially 
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from the behavior of its continuous approximation. Applying Stirling's approximation to 
(O and introducing new variables x = n\ — n%, v = p\ — P2 — let us call them "space" 
and "velocity", respectively — we obtain 

P(M),P(^^/Ie Xp {4(-^) 2 }. (3, 

This is the fundamental solution of the heat (also known as diffusion or Fokker-Planck 
equation: 

dP(x,t) | v dP(x,t) _ (i-v 2 )d 2 P(x,t) 

dt dx 2 dx 2 

Note that expression ([3]) — due to the velocity limits — 1 < v < 1 in our model — contains 

- 

wave equation 



"relativistic" fragment x . v ^ . Note also that at \v\ = 1 equation © is reduced to the 



dP(x,t) ± dP(x,t) = Q 

Now let us set a problem as is typical in mechanics: find extremal trajectories connect- 
ing two fixed points (0,0) and (X,T). As a version of the "least action principle", we 
adopt here the search of trajectories of maximum probability. The probability of trajectory 
connecting the points (0, 0) and (X, T) and passing through some intermediate point (x, t) 
is the following conditional probability 

P(x,t)P(X - x,T -t) 



P(0,0)^(x,t)^(X,T) P(X,T) 



i!(T-i)!(^)!(^)! 



The conditional probability computed for approximation (O takes the form 



(6) 



2 



T (Xt- xT) 

P m ^,t^(x,T) - V , {1 _ v 2 mT _ t) ex Pj" 2( i - v 2 )tT {T-t)} ■ (7) 

One can see essential differences between © and |7]): 

• exact probabilities © do no? depend on the velocity f (or, equivalently, on the prob- 
abilities pi, p2 of a single trial), whereas © contains artificial dependence, 

• it is easy to check that expression © allows many trajectories with the same max- 
imum probability, whereas extremals of ([7]) are deterministic trajectories, namely, 
straight lines x = ^-t. This is a typical example of emergence of deterministic be- 
haviour as a result of the law of large numbers approximation. 



7 The name of the equation depends on interpretation of the function P (x, t). 
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2.1.2. Space Symmetries in More Detail. 

A space X in our models has the structure of a graph. Graphs — we shall call them also 
lattices — are sufficient for all our purposes. In particular, they are adequate to introduce 
gauge and quantum structures. The symmetry group of the space X is the graph automor- 
phism group G = Aut (X). The automorphism group of a graph with n vertices may 
have up to n\ elements. Nevertheless, the most efficient currently algorithm designed by 
B. McKay jTOl determines the graph automorphisms by constructing compact set (no more 
than n — 1 elements, but usually much less) of generators of the group. 

Very often dynamics of a model is expressed in terms of rules defined on the neigh- 
borhoods of lattice vertices. For this sort of models with locally defined evolution rules — 
typical examples are cellular automata and the Ising model — the above mentioned group 
of local symmetries Gi oc is essential. Local rules are defined on orbits of Gi oc on edges 
from the neighborhoods of points x. Fig. [2] shows the symmetry groups G and Gi oc < G 
for some carbon and hydrocarbon molecules. 

Tetraheclrane C4//4 

G = Sym(4) 
Gioc = Do = Sym(3) 

Fullerene Cm Toric graphene n x m 

G = Z-2 x Alt(5) G = D„ x D 2m G = (Z x Z) » D 6 

G/oc = Z-i Gioc = Z-i > m ~* 00 Gi„c = D e 

n = 6 



Cubane C%H S Dodecahedrane Cai-f^o 

G = Z 2 x Sym(4) G = Z 2 x Alt(5) 






Figure 2. Symmetries of 3-valent (hydro)carbon nanostructures. 

Let us consider the role of the local group Gi oc in more detail using the buckyball f 
as an example. The incarnations of this 3-valent graph include in particular: 

• the Caley graph of the icosahedral groups Alt(5) (in mathematics); 

• the molecule of fullerene Cqq (in carbon chemistry). 



8 The classical book by F. Klein 1111 is devoted entirely to this group. 
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The symmetry group of the buckyball is G = Aut (X) = Z2 x Alt (5). The neighborhood 
of a vertex Xj takes the form . The stabilizer of Xj is Gi oc = Stabc (xj) = Z2. 

The set of neighborhood edges contains three elements: 

Ei = {e lti = (xi, x lti ) , e 2 ,j = (x*, x 2 ,i) , e 3>i = (x*, x 3ji )} . 
The set of orbits of Gi oc on consists of two orbits: 

= = {ei,i,e 2 ,i} = {e3,i}} , 

i.e., the stabilizer does not move the edge (xj,x 3j j) and swaps (x^xi^) and (xj,X2,i) . 
This asymmetry results from different roles the edges play in the structure of the buckyball: 
(xj, xij) and (xj, X2,i) are edges of a pentagon adjacent to Xj , whereas (xj, X3 5 j) separates 
two hexagons; in the carbon molecule Cqo the edge (2^,0:3^) corresponds to the double 
bond, whereas others are the single bonds. 

Naturally formulated local rules determining behavior of a system must respect decom- 
positions of neighborhoods into the orbits of the group of local symmetries. For example, 
the Hamiltonian of the Ising model on the buckyball must depend on two, generally differ- 
ent, coupling constants J\ 2 and J3. Moreover, the coupling constants may be of different 
types — feiTomagnetic or antiferromagnetic — and this may lead to interesting behavior of 
the model. Such natural Hamiltonian should take the form 

Hbucky = ~ n ^2 Si t Jl2 ( Sl >* + S2 ^ + J 3«3,j] - B ^ S i ) ( 8 ) 
i i 

where Sj, si i, s 2j i, 53,1 E X = { — 1, 1} . In a similar way the local rule for a cellular au- 
tomaton on the buckyball must have the form 

•Ej — f (^i) ^1,1) •&2,ii -^3,0 j 

where function / must be symmetric with respect to variables x\^ and X2,i, i.e., 

/ (Xj, X\ i, % 2 ,ii X3 t i) = f (Xj, X2JI, Xi j, Xs^i*) . 



2.1.3. Unification of Space and Internal Symmetries. 

Having the groups G and T acting on X and S, respectively, we can combine them into a 
single group W < Sym (£ x ) which acts on the states S = £ x of the whole system. The 
group W can be identified, as a set, with the Cartesian product T x ® G, where T x is the set 
of T-valued functions on X. That is, every element u E W can be represented in the form 
u = (a(x), a) , where a(x) E T x and a E G. A priori there are different possible ways 
to combine G and F into a single group. So selection of possible combinations should be 
guided by some natural (physical) reasons. General arguments convince that the required 
combination W should be a split extension of the group G by the group T x . In physics, it 
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is usually assumed that the space and internal symmetries are independent, i.e., W is the 
direct product T x x G with action on E x and multiplication rules: 

a(x) (a (x), a) = a(x)a(x) action, 
(a(x), a)(/3(x), b) = (a(x)/3(x), ab) multiplication. (9) 

Another standard construction is the wreath product T^xG having a structure of the semidi- 
rect product r x x G with action and multiplication 

a(x) (a (x), a) = a (ia _1 ) a (xa -1 ) , 
(a{x), a)(/3(x), b) = (a (x) /3 (xa) , ab) . (10) 

These examples are generalized by the following 

Statement. There are equivalence classes of split group extensions 

1^T X ^W^G^1 (11) 

determined by antihomomorphism^\ p, : G — > G. The equivalence is described by arbitrary 
function k : G — > G. The explicit formulas for main group operations — action on S x , 
multiplication and inversion — are 

a(x) (a (x), a) 
(a(x), a) (/3(x), b) 

(a(x), a)" 1 

This statement follows from the general description of the structure of split extensions 
of a group G by a group H: all such extensions are determined by the homomorphisms 
from G to Aut (H) (see, e.g., |[T2l ). Specializing this description to the case when H is 
the set of T-valued function on X and G acts on arguments of these functions we obtain 
our statement. The equivalence of extensions with the same antihomomorfism /i but with 
different functions k is expressed by the commutative diagram 



= a (xfi[a)) a (xK{a)) , (12) 
= (a (xK(ab)- l {i(b)K,(a)) (3 (xK(ab)~ l n(b)) , ab) , (13) 

= (a (xk (a" 1 )' 1 /-1(a)" 1 K(afj , a' 1 J . (14) 



-> w 



K 



■* w 



G 



-> G 



-> 1 



-> 1 



(15) 



where the mapping K takes the form K : (a(x), a) i-> (a (xk(o)), a) . 

Note that the standard direct and wreath products are obtained from this general con- 
struction by choosing antihomomorphisms p(a) = 1 and p(a) = a" 1 , respectively. As to 
the arbitrary function k, the choices «(a) = 1 and n(a) = a -1 , respectively, are generally 
used in the literature. 



'The term 'antihomomorphism' means that jj,(a)fi(b) = fi(ba). 
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In our computer programs (written in C) the group W is specified by two groups G and 
T and two functions fj,(a) and n(a) implemented as arrays. It is convenient in computations 
to use the following specialization: fi(a) = a~ m and = a k . For such a choice formulas 
QS-DIltake the form 

a(x) (a (x), a) 
{a(x), a) {j3 (x), b) 
(a(x), a)" 1 

Here k is arbitrary integer, but m is restricted only to two values: m = and m = 1, i.e., 
such specialization does not cover other than, respectively, direct and wreath types of split 
extentions. On the other hand, the antihomomorphisms /x(a) = 1 and fi(a) = a~ l exist 
for any group, while others depend on the particular structure of a group. Note that actions 
of G on any function f(x) are called trivial and natural for fi(a) = 1 and fi(a) = a -1 , 
respectively. 

3. Structural Analysis of Discrete Relations 

The methods of compatibility analysis, such as the Grobner basis computation or reduction 
to involutive form, are widely used to study systems of polynomial and differential equa- 
tions. In this section we develop similar techniques for discrete systems, in particular, for 
cellular automata. 

Let us consider the Cartesian product S n = Si x £ 2 e., the set of ordered 

n-tuples (o"i,o"2, . . . ,<r n ) , with oi G £j for each i. By definition, n-ary relation is any 
subset of the n-dimensional hyperparailelepiped S n . We assume that £, are finite sets of 
qi = | £j | elements that we shall call states. 

We can treat n dimensions of the hyperparailelepiped S" as elements of a set of points 
X = {x\,X2, ■ ■ ■ ,x n }. To make this initially amorphous set into a "space" (or "space- 
time") we should provide X with a structure determining how "close" to each other are 
different points. The relevant mathematical abstraction of such a structure is an abstract 
simplicial complex. The natural concept of space assumes the homogeneity of its points. 
This means that there exists a symmetry group acting transitively on X, i.e., providing 
possibility to "move" any point into any other. The homogeneity is possible only if all Sj 
are equivalent. Let us denote the equivalence class by S. We can represent £ canonically 
in the form E = {0, . . . , q — 1} , q = 

If the number of states is a prime power, q = p m , we can additionally equip the set S 
with the structure of the Galois field ¥ q . Using the functional completeness of polynomials 
— this means that any function can be represented as polynomial — over finite fields |[T3l . 
we can represent any fc-ary relation on £ as a set of zeros of some polynomial belonging to 
the ring ¥ q [x\, . . . , Xk}- Thus, the set of relations can be regarded as a system of polynomial 
equations. Although this description is not necessary (and does not work, if £j are different 
sets or q is not prime power), it is useful due to our habit to employ polynomials wherever 



a (xa m ) a y 
a lx(ab) 



k—mk+m 



(16) 

13 (x{ab)~ k b k y ab) , (17) 



a xa 



2k+m 



, a 



(18) 
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possible and capability of applying different advanced tools of polynomial algebra, such as, 
for example, the Grobner bases. 

An abstract simplicial complex (see, e.g., lfl4l ) K = (X, A) is determined by a set 
of points X = {xi,X2, ■ ■ ■ ,x n } and an assembly A of subsets of X, which are called 
simplices, such that (a) for all Xi G X {xi} G A and (b) if r C 5 G A, then r G A. The 
subsets of a simplex — they are also simplices due to (b) — are called faces. Condition 
(a) means that all one-element subsets are simplices. Clearly, the structure of the complex 
K, i.e., the set A, is uniquely determined by the simplices that are maximal by inclusion. 
Dimension of a simplex 6 is the number dim<5 = \S\ — 1. This definition is motivated 
by the fact that k + 1 points immersed in the general position into the Euclidean space 
IR n - fc form a /c-dimensional convex polyhedron. The dimension of a complex K is defined 

as the maximum dimension of all simplices in K: dim K = max dim 5. From the point 

SeA 

of view of abstract combinatorial topology, no matter how the complex can be immersed 
into the space R" — it is essential only how its simplices are connected with each other. 
However, it follows from the Nobeling-Pontryagin theorem that any (locally finite) abstract 
fc-dimensional complex can be geometrically realized in the space M 2fc+1 . We will show 
below that, for any ra-ary relation R C E n , one can regularly and uniquely construct some 
abstract simplicial complex. 

3.1. Basic Definitions and Constructions 

In addition to /c-simplices, which are singled out sets of k + 1 points, we need to consider 
arbitrary sets of point. For brevity, we shall call sets containing k points by /e-sets. Dealing 
with systems of relations defined on different sets of points, it is necessary to establish 
correspondence between the points and dimensions of the hypercube E fc . This is achieved 
by using exponential notation. The notation fixes E as the set of values of the point 
X{. For the /c-set S = {x\, . . . , x^}, we introduce the notation E* 5 = E^' 1 ^ x • • • x E^*^. 
The set 5 is called the domain of the relation R s . We will call the whole hypercube E 5 a 
trivial relation. Accordingly, R s C E 5 denotes a relation given on the set of points 5. 

3.1.1. Relations. 

Thus, we have: 

Definition 1 (relation). A relation R s on the set of points 6 = {x\, . . . ,x^\ is any 
subset of the hypercube E 5 ; i.e., R 5 C T, 5 . 

The relation R s can be regarded as the Boolean- valued function R s : E 5 — > {0, 1}. We can 
think of variables taking values in E and write the relation as 

a = R s (xt, ...,x k ), a G {0, 1} . 

An important special case of relations: 

Definition 2 (functional relation). A relation R s on the set of points 
5 = {xi, . . . , Xk} is called functional if there is a position i G (1, . . . , k) such that 
for any a±, . . . , fTj_i, cr i+ i, . . . ,a k ,q,T G E 

from (<ti, . . . ,o-i-i,s,o- i+1 , . . . ,0-fc) G R 5 and (01, . . . , cr^-i, r, a i+ i, . . . , a k ) G R s 
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it follows that ? = r. 

In terms of variables the functional relation R s can be written in the form 

Xi = F(x u .. . ,Xi_i,x i+ i,... ,x k ) , where F : T, 8 \ {xi} -> S. 

We need to be able to extend relations from subsets of points to larger sets: 
Definition 3 (extension of relation). For given set of points 5, its subset r C 5 and 
relation R T on the subset r, we define the extension of R T as the relation 

R s = R T x S 5 \ T . 

This definition, in particular, allows the relations R Sl , . . . , R Sm defined on different do- 
mains to be extended to the common domain, i.e., to the union Si U ■ ■ ■ U 5 m . 
Logical implications of the relations are defined in a natural way: 
Definition 4 (consequence of relation). A relation Q s is called a consequence of the 
relation R 5 if R 5 C Q 5 C T, 5 ; i.e., Q 5 is arbitrary superset of the set R 5 . 

The relation i? 5 may have many different consequences: their total number (including 

r s I y><5 I I d<5 I 

it itself and the trivial relation £ ) is evidently equal to 21 I I I. 

It is natural to single out the consequences that can be reduced to relations on smaller 
sets of points: 

Definition 5 (proper consequence). A nontrivial relation Q T is called the proper conse- 
quence of the relation R s if t is a proper subset of 5 (i.e., r C 6) and the relation Q T x X 5 \ T 
is a consequence of R s . 

We call relations that have no proper consequences the prime relations. 

3.1.2. Compatibility of Systems of Relations. 

The compatibility of a system of relations can naturally be defined by the intersection of 
their extensions to the common domain: 

Definition 6 (base relation). The base relation of the system of relations 
R Sl , . . . , R 5m is the relation 

m m 

R s = p| R 5 > x E Si \ s , where S = (J 
i=i i=i 

Let us make two comments for the polynomial case q = p n , where the standard tool for the 
compatibility analysis is the Grobner basis method: 

• The compatibility condition determined by the base relation can be represented by a 
single polynomial, unlike the Grobner basis, which is normally a system of polyno- 
mials. 

• Any possible Grobner basis of polynomials representing the relations R Sl , . . . , R Sm 
corresponds to some combination of consequences of the base relation. 
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3.1.3. Decomposition of Relations. 

If a relation has proper consequences, we can try to express it as far as possible in terms of 
these consequences, i.e., relations on smaller sets of points. To this end we introduce 

Definition 7 (canonical decomposition). The canonical decomposition of a relation R s 
with proper consequences Q Sl , . . . , Q Sm is the relation 

*w n (n ® 5i x e<5V5< j ' (i9) 

where the factor P s is defined by the following 

Definition 8 (principal factor). The principal factor of the relation R s with proper 
consequences Q Sl , . . . , Q Sm is the relation 



P 5 = R s (J It, 5 \ p| Q 5 * x . 



\ i=l / 

The principal factor is the maximally "free" — i.e., the closest to the trivial — relation that, 
together with the proper consequences, makes it possible to recover the initial relation. 

If the principal factor in the canonical decomposition is trivial, the relation is completely 
reduced to relations on smaller sets of points. 

Definition 9 (reducible relation). A relation R s is said to be reducible if it can be 
represented as 

m 

R s = f] Q Si x (20) 

i=l 

where 5- L are proper subsets of S. 

This definition makes it possible to impose a "topology" — i.e., the structure of an 
abstract simplicial complex with the corresponding theories of homologies, cohomologies, 
etc. — on an arbitrary n-ary relation R C S n . This is achieved by 

• naming the dimensions of the hypercube S" as the "points" x\, . . . , x n € X, 

• decomposing R (which can now be denoted by R x ) into irreducible components, 

• and defining the maximal simplices of the set A as the domains of irreducible com- 
ponents of the relation R x . 



3.1.4. On Representation of Relations in Computer. 

A few words are needed about computer implementation of relations. To specify a k-ary 
relation R k we should mark its points within the /c-dimensional hypercube (or hyperparal- 
lelepiped) T, k , i.e., define a characteristic function \ '■ ^ k ~ * {0> 1} , with x{&) = 1 or 
according as a € R k or a ^ R k . Here a = (o"o,o"i, . . . ,<7fc_i) is a point of the hy- 
percube. The simplest way to implement the characteristic function is to enumerate all the 
q k hypercube points in some standard, e.g., lexicographic order. Then the relation can be 
represented by a string of q k bits ao«i • • • a ? fe_i in accordance with the table: 
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We call this string bit table of relation. Symbolically Bit Table [ig\ := (<r G R k ) . Note 
that a is the ("little-endian") representation of the number ig in the radix q: 

iff = 00 + o\q H h <7»g l H h o-fc^i^'" 1 . 

In the case of hyperparallelepiped X fc = Xi x £2 x • • • x Sfc one should use the multi-radix 
representation of integers: 

iff = 0"o + 01 X <?lH 1- fi x 9i 92 ■ • • 9i H h o-fc-i x g x g 2 • ■ ■ qk-l, 

where < < ft+i, i £ [0, . . . , k — 1], 

The characteristic function (bit table) can be represented as the binary integer 

x = ao + ai2 + • • • + + • • • + a^^ 9 * -1 . (21) 

Most manipulations with relations ai - e reduced to very efficient bitwise computer com- 
mands. Of course, symmetric or sparse (or, vice versa, dense) relations can be represented 
in a more economical way, but these are technical details of implementation. 

3.2. Illustration: Application to Some Cellular Automata 
3.2.1. J. Conway's Game of Life. 

The "Life family" is a set of 2-dimensional, binary (i.e., S = {0, 1}; q = 2) cellular au- 
tomata similar to Conway's Life, which rule is defined on 9-cell (3x3) Moore neighborhood 
and is described as follows. A cell is "born" if it has exactly 3 "alive" neighbors, "survives" 
if it has 2 or 3 such neighbors, and "dies" otherwise. This rule is symbolized in terms of 
the "birth'V'survival" lists as B3/S23. Another examples of automata from this family are 
HighLife (the rule B36/S23), and Day&Night (the rule B3678/S34678). The site 03 con- 
tains collection of more than twenty rules from the Life family with Java applet to run these 
rules and descriptions of their behavior. 

Generalizing this type of local rules, we define a k-valent Life rule as a binary rule on a 
fc-valent neighborhood (we adopt that x\, . . . ,Xk sue neighbors of Xk+i of the central cell 
Xk+i) described by two arbitrary subsets of the set {0, 1, ... , k}. These subsets B,S C 
{0,1, ... ,k} contain conditions for the one-time-step transitions x^+i — > x' k+l of the forms 
0—>l and 1—^1, respectively. Since the number of subsets of any finite set A is 2^1 and 
different pairs BIS define different rules, the number of different rules defined by two sets 
B and S is equal to 2 k+1 x 2 k+l . Thus, the total number of fc-valent rules described by the 
"birth'V'survival" lists is 

N B/S)k = 2 2k+2 . (22) 
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There is another way to characterize this type of local rules. Let us consider /c-valent 
rules symmetric with respect to the group Sym (k) of all permutations of k outer points of 
the neighborhood. We shall call such rules fc-symmetric. It is not difficult to count the total 
number of different q-aiy ^-symmetric rules: 

N q , Sym (k) =q [ «-' > Q . (23) 

We see that (l23l evaluated at q = 2 coincides with (1221 . i.e., _/V 2) s ym (fc) = ^B/s, k- Since 
/c-valent Life rules are obviously fc-symmetric we have the following 

Proposition. For any k the set of k- symmetric binary rules coincides with the set ofk-valent 
Life rules. 

This proposition implies in particular that one can always express any /c-symmetric 
binary rule in terms of the "birth'V'survival" lists. 

The local relation of Conway's Life automaton Rq L is defined on the 10-set 5 = 

{x 1 ,...,x 10 }: 




Here the point x\q = x' g is the next-time-step of the point xg. By construction, elements of 
the 10-dimensional hypercube belong to the relation of Conway's Life automaton, i.e., 
(xi, . . . , xiq) e Rql' i n tne following cases: 

1- (ELi x i = 3) A (xiq = 1), 

2 - (Ya=1 x i = 2 ) A ( x 9 = x w), 

3. xio = 0, if none of the above conditions holds. 

The number of elements of the relation Rq L is l-R^il = ^ e relation R S CL , as is 
the case for any cellular automaton, is functional: the state of xio is uniquely determined 
by the states of other points. The state set S = {0, 1} can be additionally endowed with 
the structure of the field F2 . We accompany the below analysis of the structure of R 5 CL by 
description in terms of polynomials from F2 [xi, ■ ■ ■ , xio] • This is done only for illustrative 
purposes and for comparison with the Grobner basis method. In fact, we transform the 
relations into polynomials only for output. Transformation of any relation into polynomial 
form can be performed by computationally very cheap multivariate version of the Lagrange 
interpolation. In the case q = 2, the polynomial which set of zeros corresponds to a relation 
is constructed uniquely. If q = p n > 2, there is a freedom in the choice of nonzero values 
of constructed polynomial, and the same relation can be represented by many polynomials. 
The polynomial representing R S CL takes the form 



Pcl = x w + x 9 (n 7 + n 6 + n 3 + n 2 ) + n 7 + n 3 , 



(24) 
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where IT^ = ITfc (xi, . . . , x%) is the kth elementary symmetric polynomial defined for n 
variables x±, . . . , x n by the formula: 

life • • • ) = ^ ^ ' ' ' 

l<ii<i2<--<ifc<n 

Hereafter, we will use the following notation: 

n fc = n fc (xi, ...,x 8 ), rife = n fc (xi, ,x 8 ) , 
n fc J = n fc (xi, . . . ,fj, . . . ,£j, . . . ,x 8 ) . 

Applying the computer program to Rq L , we find that the relation R & CL is reducible and has 
the decomposition 

^=< w n(n^ K} )> (25) 

where (ii, . . . ,27) is arbitrary 7-element subset of the set (1, . . . ,8). For brevity, we 
dropped in (T25T ) the trivial factors entering into the general formula (|20b . 

The eight relations R^ Xl ^ (1 < i < 8; for decomposition (l25T ). it suffices to take any 
seven of them) have the following polynomial form: 

x 9 x 10 (ir 6 + n| + rr 2 + ni) + x 10 (rr 6 + n| + 1) + x 9 (n 7 + rr 6 + rr 3 + n 2 ) = o. 

Accordingly, the relation R^ xg ^ has the form 

xio (n 7 + n 6 + n 3 + n 2 + 1) + n 7 + n 3 = 0. 

The relations R^ x ^ and R^ are irreducible but «of prime, and can be expanded in 
accordance with formula (fl9l ). Continuing the decomposition iterations, we finally obtain 
the following system of relations (in the polynomial form) that are satisfied for Conway's 
Life: 

x 9 xi (n| + ni) +xi (n| + i) + x 9 (n 7 + n| + n| + n 2 ) = o, (26) 

xio(n 3 + n 2 + i) + n 7 + n 3 = o, (27) 

(x 9 xio + xio) (l$ + Y?i + U\ j + l) = 0, (28) 

xi (n| + n| + ni + i) = o, (29) 

Xl0^iiXi2Xi 3 Xi 4 — 0. (30) 

One can easily interpret the simplest relations (l30l : if the point xio is in the state 1, then at 
least one point in any set of four points surrounding x 9 must be in the state 0. 

The above analysis of the relation R S CL takes < 1 sec on a 1.8GHz AMD Athlon note- 
book with 960Mb. 

To compute the Grobner basis we must add to polynomial (l24l ) ten polynomials 

xf + Xi, i = l,..., 10 

corresponding to the relation x q = x that holds for all elements of any finite field ¥ q . 

Computation of the Grobner basis over F 2 with the help of Maple 9 gives the following: 
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• pure lexicographic order with variable ordering xio >- Xg >-•••>- x\ does not 
provide any new information leaving initial polynomial (l24l) unchanged; 

• pure lexicographic order with variable ordering x\ >~ X2 >- ■ ■ ■ >- xio reproduces 
relations (l26l) — (l30l (modulo several polynomial reductions violating the symmetry 
of polynomials); the computation takes 1 h 22 min; 

• degree-reverse-lexicographic order also reproduces system (l26l > — (1301) (same com- 
ment as above); the times are: 51 min for the variable ordering x\ >- x% >- ■ ■ ■ >- x\$, 
and 33 min for the ordering x\q >- xg >~ . . . y x\. 

3.2.2. Elementary Cellular Automata. 

Simplest binary, nearest-neighbor, 1 -dimensional cellular automata were named elementary 
cellular automata by S. Wolfram, who has extensively studied their properties |[T6l . A 
large collection of results concerning these automata is presented in Wolfram's online atlas 
iMTl . In the exposition below we use Wolfram's notations and terminology. The elementary 
cellular automata are simpler than Conway's Life, and we may pay more attention to the 
topological aspects of our approach. 

Local rules of the elementary cellular automata are defined on the 4-set 5 = {p, q, r, s} 

@ © © 

which can be pictured by the icon y/ • A local rule is a binary function of the form 

s = f(p, q, r). There are totally 2 23 = 256 such functions, each of which can be indexed 
with an 8 -bit binary number. 

Our computation with relations representing the local rules shows that the total number 
256 of them is divided into 118 reducible and 138 irreducible relations. Only two of the 
irreducible relations appeared to be prime, namely, the rules 105 and 150 in Wolfram's nu- 
meration. This numeration is based on the "big-endian" — i.e., opposite to our convention 
(|2TI) — representation of binary numbers. Note, that the prime rules 105 and 150 have 
linear polynomial forms: s=p + q + r + l and s = p + q + r, respectively. 

We consider the elementary automata on a space-time lattice with integer coordinates 
(x, t), i.e., x G Z or x G Z m (spatial m-periodicity), t G Z. We denote a state of the point 
on the lattice by u(x,t) G S = {0, 1}. Generally the points are connected as is shown in 
the picture 



•< X »- 




The absence of horizontal ties expresses the independence of "space-like" points in cellular 
automata. 

Reducible Automata. The analysis shows that some automata with reducible local rela- 
tions can be represented as unions of automata defined on disconnected subcomplexes: 
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Two automata and 255 are determined by unary relations s = and s = 1 on the 
disconnected set of points: 



Note that unary relations are usually called properties. 

Six automata 15, 51, 85, 170, 204 and 240 are, in fact, disjoint collections of spacially 
zero-dimensional automata, i.e., single cells evolving in time. As an example, let us 

® ® p 

consider the automaton 15. The local relation is defined on the set m ^ * ts 

bit table is 0101010110101010. This relation is reduced to the relation on the face 

® 

@ with bit table 0110. The spacetime lattice is split in the following way: 




The bit table 01 10 means that the points p and s can be only in opposite states, and 
we can write immediately the general solution for the automaton 15: 

u(x, t) = a(x — t) +t mod 2, 

where u(x, 0) = a(x) is an arbitrary initial condition. 

Each of the ten automata 5, 10, 80, 90, 95, 160, 165, 175, 245, 250 is decomposed into 
two identical automata. As an example let us consider the rule 90. This automaton 
is distinguished as producing the fractal of topological dimension 1 and Hausdorff 
dimension In 3/ In 2 « 1.58 known as the Sierpinski sieve (or gasket or triangle). Its 

© © © 

local relation on the set is represented by the bit table 101001010101 1010. 

The relation is reduced to the relation with the bit table 



10010110 on the face 



V 6 



(31) 



It can be seen from the structure of face (|3TT > that the spacetime lattice is split into 
two identical independent complexes as is shown 





U 
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To find a general solution of the automaton 90 it is convenient to use the polynomial 
form of relation (I3TT) s + p + r = 0. With this linear expression, the general solution 
is easily constructed: 

u(x,t) = '^2 (\a(x — t + 2k) mod 2, u(x,0) = a(x). 
k=o ^ ' 

Using Proper Consequences. Proper consequences — even if they are not functional — 
can provide useful information on the behavior of a cellular automaton. 

For example, 64 automata^ (with both reducible and irreducible local relations) have 
proper consequences with the bit table 

1101 (32) 

on, at least, one of the faces 




The algebraic forms of relation (l32l on faces (l33l are ps + s = 0, qs + s = 0, rs + s = 0, 
respectively. 

Relation (l32l is not functional, and hence can not describe any deterministic evolution. 
Nevertheless, it imposes severe restrictions on the behavior of the automata having such 
proper consequences. The features of the behavior resulting from relation (l32l > are clearly 
seen from many of computational results presented in the atlas iTTTll - A typical pattern 
from this atlas is reproduced in Fig. [3l where several evolutions of the automaton 168 
are presented. In the figure, 0's and l's are denoted by the empty and filled square cells, 
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Figure 3. Rule 168. Several random initial conditions 



respectively. Note that the authors of the figure use a spatially periodic condition: x G Z30. 



The complete list of these automata in Wolfram's numeration is as follows: 2, 4, 8, 10, 16, 32, 34, 40, 42, 
48, 64, 72, 76, 80, 96, 112, 128, 130, 132, 136, 138, 140, 144, 160, 162, 168, 171, 174-176, 186, 187, 190-192, 
196, 200, 205, 206, 208, 220, 222-224, 234-239, 241-254. 
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The local relation of the automaton 168 — its polynomial form is pqr + qr +pr + s = 
— has the proper consequence rs + s = 0. Relation (l32l) means that if, say r, as for the 
rule 168, is in the state 1 then s may be in both states or 1, but if the state of r is 0, then 
the state of s must be 0: 

r = l^s = 0Vs = l, 
r = =>■ s = 0. 

One can see that all evolutions in Fig. [3] consist of diagonals x = xq — t directed leftward 
and downward. Each diagonal begins with a several units, but after the first appearance of 
zero all subsequent points along the diagonal are zeros. 

Canonical Decomposition vs. Grobner Basis. In this paragraph we compare our canon- 
ical decomposition ( fl9l ) with Grobner basis in the polynomial case. Let us begin with two 
examples of elementary cellular automata. The Grobner bases are computed in the total 
degree and reverse lexicographical order of monomials. The trivial polynomials p 2 + p, 
q 2 + q, r 2 + r and s 2 + s are omitted in the Grobner bases descriptions. 

• Automaton 30 is remarkable by its chaotic behavior and is even used as a random 
number generator in Mathematica. 

Relation: 1001010101101010 or qr + s + r + q + p = 0. 
Canonical Decomposition: 

Proper consequences: 



face 




bit table 11011110 11011110 

polynomial qs + pq + q rs + pr + r. 

Principal factor: 1011111101111111 or qrs +pqr + rs + qs +pr +pq + s +p = 0. 
Grobner basis: {qr + s + r + q + p, qs + pq + q, rs + pr + r} . 

Thus for the rule 30 the polynomials of the canonical decomposition coincide (mod- 
ulo obvious polynomial substitutions) with the Grobner basis. 

• Automaton 110 is, like a Turing machine, universal, i.e., it can simulate any compu- 
tational process, in particular, any other cellular automaton. 

Relation: 1100000100111110 or pqr + qr + s + r + q = 0. 

Canonical Decomposition: 

Proper consequences: 
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bit table 11011111 11011111 10010111 

polynomial pqs + qs + pq + q prs + rs + pr + r qrs + s + r + q. 

Principal factor: 1111111111111110 or pqrs = 0. 
Grobner basis: 

{prs + rs + pr + r, gs + rs + r + q, qr + rs + s + q, pr + pq + ps} . 

For automaton 1 10, the polynomials of the Grobner basis are not identical with the 
polynomials of the canonical decomposition. The system of relations defined by the 
Grobner basis is: 

R {p,r,s} = 1101im = ( prs + rs + pr + r = ), 

Rl q ' r ' s} = 10011111 = (qs + rs + r + q = 0), 

R {<l,r,s} = 10110111 = (qr + rs + s + q = 0), 

R {p,q,r;s} = moioillOlllllO = {pr + pq + ps = 0) . 

In general, the folowing differences between our approach and the Grobner basis method 
can be mentioned. 

• In contrast to a Grobner basis, a base relation, defined as intersection of conditions, 
agrees with the standard in logic and set theory notion of compatibility. 

• In contrast to a canonical decomposition a Grobner basis may look beyond the poly- 
nomial context as a collection of accidental supersets. 

• There is some analogy between Grobner bases and canonical decompositions — in 
fact, they coincide in about half of cases in our computations. 

• Canonical decomposition is more efficient for problems with polynomials of arbitrary 
degree — the above computation with Conway's automaton is an example. 

• For small degree problems with large number n of indeterminates the Grobner basis 
outperforms canonical decomposition — the number of polynomials of bounded de- 
gree is a polynomial function of n, whereas the algorithm of canonical decomposition 
scans exponential number q n of the hypercube points. 

4. Soliton-like Structures in Deterministic Dynamics 

Symmetries of deterministic systems impose severe restrictions on the system dynamics 
loll . In particular, for the first ordeiP^I functional relations: 

"This means that evolution relation (fJJ takes the form St = F (s t _i). 



Analysis of Discrete Dynamical Systems 



23 



• dynamical trajectories pass group orbits in non-decreasing order of orbit sizes, 

• periodic trajectories lie within orbits of the same size. 

One of the characteristic features of dynamical systems with non-trivial symmetries is for- 
mation of moving form-preserving structures. 

Let us begin with a simple example. Consider a cube X whose vertices take values 
in two-element set, say £ = {0, 1} . By the way, as is clear from Fig. 0] a cube can be 
interpreted as a simplest "finite model of graphene". The 48-element symmetry group of a 



a a b a. 




b aba 



Figure 4. The graph of cube forms 4-gonal (6 tetragons) lattice in sphere § 2 and 6-gonal (4 
hexagons) lattice in torus T 2 . 



cube has the structure G = Z2 x Sym (4) . The group is generated by 3 elements: 



1. 120° rotation around diagonal of the cube; 



2. 90° rotation around axis passing through the centers of opposite cube faces; 

3. reflection interchanging opposite faces of the cube. 

Total number of states of the model is |S X | = 2 s = 256. If we assume that the group V is 
trivial, then W = T x xG = lxG = G. The group W splits the set S x into 22 orbits in 



accordance with the table: 



Size of orbits 



1 2 



12 24 



2 12 2 5 



Number of orbits 

Let us consider a deterministic dynamical system on the cube, namely, symmetric bi- 
nary 3-valent cellular automaton with the rule 86. The number 86 is the "little endian" 
representation of the bit string 01101010 taken from the last column of the rule table with 
Sym (3)-symmetric combinations of values for x± i i,X2,i, xsi 
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Here Xi is value of ith vertex of the cube; Xi,j, X2,ijX3 ) i are values of the cube vertices 
adjacent to the ith one and x\ is the next time value of ith vertex. The rule can also be rep- 
resented in Conway's Life style "Birth7"Survival" notation as B123/S0, or as polynomial 
over the field F2 

4 = Xi + Il3 + n2 + ni, 

where 111 = x lti + x 2 ,i + x 3ti , U 2 = x lti x 2 ,i + x lti x 3 ^ + X2,jX 3) i, 1I3 = x lti x 2 ,iX 3ti are 
elementary symmetric functions. 

The phase portrait of the automaton is shown in Fig. [5j where the group orbits are 
represented by circles containing the ordinal number^] of orbits within. The numbers over 
orbits and within cycles are sizes of the orbits (recall that all orbits belonging to the same 
cycle have equal sizes — see the beginning of this section). The rational number p indicates 
the weight of the corresponding element of the phase portrait. In fact, p is a probability for 
randomly chosen state to appear in an isolated cycle or to be caught by an attractor: p = 
(size of basin)/ {total number of states). Here size of basin is sum of sizes of orbits involved 
in the struture. 




Figure 5. Rule 86. Equivalence classes of trajectories on hexahedron. 



Generalizing this example, we see that if the symmetry group W splits the state set 
S x of deterministic dynamical system into finite number of orbits, then after some lapse of 
time any trajectory comes inevitably to a cycle over some finite sequence of orbits. This 
just means formation of soliton-like structures. Namely, let us consider evolution 

a to (x) -)• a tl (x) = A tlto (ot (x)) . (34) 

If the states at the moments io an d t\ belong to the same orbit: o- to {x) G Oi and o~t {x) G 
Oi, Oi C S x ; then evolution (l34l can be replaced by the group action 

a tl (x) = a to (x)w, w G W, 



12 These numbers are specified by the computer program in the course of computation. 
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i.e., the initial state ("shape") at (x) is reproduced after some "movement" in the space £ . 

The following are several examples (including continuous cases) of cycles over group 
orbits: 

• traveling waves a(x — vt) in mathematical physics — the Galilei group; 

• "generalized coherent states" in quantum physics — unitary representations of com- 
pact Lie groups; 

• "spaceships" in cellular automata — lattice symmetries. 

Let us consider the "glider" — one of the "spaceships" in Conway's Life automaton. 

The space X of Conway's Life is a square lattice. For the finiteness, we shall assume 
that the lattice is closed into the N x N torus. In the general case N ^ 4 the symmetry 
group of X is the semidirect product of two-dimensional translations T 2 = Zjy x Z^r and 
the dihedral group Ds = Z4 x Z 2 : 

G = T 2 x D 8 , if N = 3,5,6, ... ,00. (35) 

In the case N = 4 the translation subgroup T 2 = Z4 x Z4 is not normal and G has a bit 
more complicated structure Q : 

normal closure of T 2 

G = ((((Z 2 x D 8 ) x Z 2 ) x Z 3 ) x Z 2 ) xZ 2 . (36) 

The extra symmetry Z3 in (l36l can be explained by the Z3 symmetry of the four-vertex 



Dynkin diagram D4 = 




associated with the case iV = 4. 



The set of local (cell) states of Conway's Life is £ = {"dead", "alive"} = {0, 1}. Since 
the local rule of Conway's Life is not symmetric with respect to the transposition 0f> lof 
the local states, the internal symmetry group is trivial, i.e., Y = {1} and hence T x = {1}. 
Thus, we have W = T x xG = lxG = G. The natural action of W on functions 
a{x) G £ x takes the form o(x)w = a (xg~ l ), where w = (1, g), g € G. 

Fig. [6] shows four steps of evolution of the glider. The figure demonstrates how the 
evolution is reduced to the group action. N > 4 is assumed. 

£71 (T 2 CJ ?, °4 °5 
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Figure 6. Example of soliton-like structure. "Glider" in Conway's Life is cycle in two 
orbits of the group G = T 2 x D§: configurations as and a 4 are obtained from o~\ and 
cr 2 , respectively, by the same combination of downward shift, 90° clockwise rotation and 
reflection in respect to vertical. 
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Comments on Reversibility in Discrete Systems. 

A typical deterministic dynamical system is irreversible — it's phase portrait modulo group 
orbits looks like in Fig. [5] We see there several isolated and limit cycles (fixed points are 
regarded as cycles of unit length) accompanied by influxes flowing into the limit cycles. In 
contrast to continuous systems, any discrete system "forgets" influxes after some time and 
appears in either isolated or limit cycles. After loss of information about influxes both types 
of cycles became physically indistinguishable and the system behaves just like reversible. 
This might be a hint for explanation of observable reversibility of the fundamental laws of 
nature. 

In this connection we would like to mention recent works of G. 't Hooft. One of the 
difficulties of quantum gravity is a conflict between irreversibility of gravity — information 
loss at the black hole horizon — with reversibility and unitarity of the standard quantum 
mechanics. In several papers of recent years (see, e.g., IPT81 [191 ) 't Hooft developed an 
approach to reconciling both theories. The approach is based on the following assumptions 

• physical systems have discrete degrees of freedom at tiny (Planck) distance scales; 

• the states of these degrees of freedom form primordial basis of Hilbert space (with 
nonunitary evolution); 

• primordial states form equivalence classes: two states are equivalent if they evolve 
into the same state after some lapse of time; 

• the equivalence classes by construction form basis of Hilbert space with unitary evo- 
lution described by time-reversible Schrodinger equation. 

In our terminology this corresponds to transition to limit cycles: in a finite time of evolution 
the limit cycle becomes physically indistinguishable from reversible isolated cycle — the 
system "forgets" its pre-cycle history. 

This type of irreversibility hardly can be observed experimentally (assuming, of course, 
that considered models may have at all any relation to physical reality). The system should 
probably spend time of order the Planck unit (?» 10 -44 sec) out of a cycle and potentially 
infinite time on the cycle. Nowadays, the shortest experimentally fixed time is about 10 -18 
sec or 10 26 Planck units. 

5. Mesoscopic Lattice Models 

Discrete symmetry analysis simplifies manipulations with microcanonical ensembles and 
search of phase transitions. This allows to reveal subtle details in behavior of mesoscopic 
models. 

5.1. Statistical Mechanics 

As we mentioned earlier, the state of deterministic dynamical system at any point of time is 
determined uniquely by previous states of the system. A Markov chain — for which tran- 
sition from any state to any other is possible with some probability — is a typical example 
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of non-deterministic dynamical system. In this section we apply symmetry approach to the 
lattice models in statistical mechanics. These models can be regarded as special instances 
of Markov chains. Stationary distributions of the Markov chains are studied by the methods 
of statistical mechanics. 

The main tool of conventional statistical mechanics is the Gibbs canonical ensemble 
— imaginary collection of identical systems placed in a huge thermostat with temperature 
T. The statistical properties of canonical ensemble are encoded in the canonical partition 
function 

Z= £ e ~ E ^ ksT . (37) 

Here S x is the set of microstates, E a is energy of microstate a, k E is Boltzmann's con- 
stant. The canonical ensemble is essentially asymptotic concept: its formulation is based 
on approximation called "thermodynamic limit". For this reason, the canonical ensemble 
approach is applicable only to large (strictly speaking, infinite) homogeneous systems. 

5.2. Mesoscopy 

Nowadays much attention is paid to study systems which are too large for a detailed mi- 
croscopic description but too small for essential features of their behavior to be expressed 
in terms of classical thermodynamics. This discipline — often called mesoscopy — covers 
wide range of applications from nuclei, atomic clusters and nanotechnological structures to 
multi-star systems ll20ll2Tll22l . To study mesoscopic systems one should use more funda- 
mental microcanonical ensemble instead of canonical one. A microcanonical ensemble is a 
collection of identical isolated systems at fixed energy. Its definition does not include any 
approximating assumptions. In fact, the only key assumption of a microcanonical ensemble 
is that all its microstates are equally probable. This leads to the entropy formula 

S E = k B In n E , (38) 

or, equivalently, to the microcanonical partition function 

Q E = e SE/kB . (39) 

Here Q E is the number of microstates at fixed energy E : Y, M E = |S X |. In what follows 

E 

we will omit Boltzmann's constant assuming k B = 1. Note that in the thermodynamic limit 
the microcanonical and canonical descriptions are equivalent and the link between them is 
provided by the Laplace transform. On the other hand, mesoscopic systems demonstrate 
experimentally and computationally observable peculiarities of behavior like heat flows 
from cold to hot, negative specific heat or "convex intruders" in the entropy versus energy 
diagram, etc. These anomalous — from the point of view of canonical thermostatistics — 
features have natural explanation within microcanonical statistical mechanics ll22l . 

5.2.1. Lattice Models. 

In this section we apply symmetry analysis to study mesoscopic lattice models. Our ap- 
proach is based on exact enumeration of group orbits of microstates. Since statistical studies 
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are based essentially on different simplifying assumptions, it is important to control these 
assumptions by exact computation, wherever possible. Moreover, we might hope to reveal 
subtle details in behavior of system under consideration with the help of exact computation. 

As an example, let us consider the Ising model. The model consists of spins placed on 
a lattice. The set of vertex values is £ = {— 1,1} and the interaction Hamiltonian is given 

by 

H = -jJ2siSj-Bj2si, (40) 

(m) » 

where Si,Sj £ S; J is a coupling constant (J > and J < correspond to ferromag- 
netic and antiferromagnetic cases, respectively); the first sum runs over all edges of 
the lattice; B is an external "magnetic" field. The second sum M = Yl s i * s called the 

i 

magnetization. To avoid unnecessary technical details we will consider only the case J > 
(assuming J = 1) and B = in what follows. 

Let us remind that if the local symmetry group Gi oc decomposes the sets of edges 
of lattice neighborhoods into nontrivial orbits, then the interaction Hamiltonian should be 
modified (see, e.g., Eq. © on page©. 

Since Hamiltonian and magnetization are constants on the group orbits, we can count 
numbers of microstates corresponding to particular values of these functions - and hence 
compute all needed statistical characteristics - simply by summation of sizes of appropriate 
orbits. 

Fig. |7]shows microcanonical partition function for the Ising model on the dodecahedron 



Here total number of microstates |S X | = 1048576, number of lattice vertices Nx = 20, 
energy E is value of Hamiltonian. 

Of course, other characteristics of the system can be computed easily in this way. 

5.3. Phase Transitions 

Needs of nanotechnological science and nuclear physics attract special attention to phase 
transitions in finite systems. Unfortunately classical thermodynamics and the rigorous the- 
ory of critical phenomena in homogeneous infinite systems fails at the mesoscopic level. 
Several approaches have been proposed to identify phase transitions in mesoscopic sys- 
tems. Most accepted of them is search of "convex intruders" ll23l in the entropy versus 
energy diagram. In the standard thermodynamics there is a relation 



d 2 s 



dE 2 



where Cy is the specific heat at constant volume. 

Relation (|4TI ) implies that d 2 S/dE 2 \ v < and hence the entropy versus energy di- 
agram must be concave. Nevertheless, in mesoscopic systems there might be intervals of 
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0.15 




Figure 7. Ising model on dodecahedron. Microcanonical distribution. 



energy where d 2 S/dE 2 \ > 0. These intervals correspond to first-order phase transitions 
and are called "convex intruders". From the point of view of standard thermodynamics 
one can say about phenomenon of negative heat capacity, of course, if one accepts that it 
makes sense to define the variables T and Cy as temperature and the specific heat at these 
circumstances. In ll24ll it was demonstrated via computation with exactly solvable lattice 
models that the convex intruders flatten and disappear in the models with local interactions 
as the lattice size grows, while in the case of long-range interaction these peculiarities sur- 
vive even in the limit of an infinite system (both finite and long-range interacting infinite 
systems are typical cases of systems called nonextensive in statistical mechanics). 
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Figure 8. Ising model on dodecahedron. "Convex intruders" on entropy vs. energy diagram 
indicate mesoscopic phase transitions. 
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A convex intruder can be found easily by computer for the discrete systems we discuss 
here. Let us consider three adjacent values of energy I^-i, Ei, and corresponding 
numbers of microstates fi/r , fi a , fip.., • In our discrete case the ratio El+1 ~ El is al- 
ways rational number p/g and we can write the convexity condition for entropy in terms of 
numbers of microstates as easily evaluated inequality 

n F +q <n p F n q E . (42) 

As a rule E i+ i — E^ = E^ — and inequality (1421 takes the form 

This form means that within convex intruder the number of states with the energy Ei is less 
than geometric mean of numbers of states at the neighboring energy levels. 

Fig. [8] shows the entropy vs. energy diagram for the Ising model on dodecahedron. The 
diagram has apparent convex intruder A in the specific energy interval [—1.2,-0.9]. Exact 
computation reveals also a subtle convex intruder B in the interval [—0.8,-0.6]. 



6. Gauge Connection and Quantization 

All most successful contemporary theories in fundamental physics are gauge theories. 
There are also numerous applications of gauge theories in mathematics (topological quan- 
tum field theory, invariants of 3- and 4-manifolds, monoidal categories, Hopf algebras and 
quantum groups, etc. IT25T0 . 

In fact, the gauge principle expresses the very general idea that in spite of the fact that 
any observable data are represented in different "reference frames" at different points in 
spaced and time, there should be some way to compare these data. 

6.1. Discrete Gauge Principle 

At the set-theoretic level, i.e., in the form suitable for both discrete and continuous cases, 
the main concepts of the gauge principle can be reduced to the following. We have 

• a set T, discrete or continuous time, T = Z or T = M; 

• a set X, space; 

• the sets T and X are combined into a space-time M. = X x T; 

• a set E, local states; 

• a group T < Sym (£) acting on X, internal symmetries; 

• identification of data describing the states from £ makes sense only modulo symme- 
tries from T — this is arbitrariness in the choice of a "reference frame"; 



~ Consideration only time evolution of general set of states <S leads to the trivial gauge structures. Gauge 
theories of interest are possible if there exists underlying space structure, i.e., S = E x . 
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• there is no a priori connection between data (i.e., between reference frames) at dif- 
ferent points x, y G M. — we should impose this connection (or parallel transport) 
explicitly as T-valued function on edges (pairs of points) of abstract graph: 

q(y) = a(x)n (x, y) , ir(x,y)er, a(x),^(y) G S; 

the connection ir (x, y) has the obvious property tt (x, y) = ir (y, x)" 1 ; 

• a connection tt (x, y) is called trivial if it can be expressed in terms of a function on 

vertices of the graph: tt (x, y) = p (x) p , p (x) ,p(y) G T; 

• invariance with respect to the gauge symmetries depending on time and space leads 
to the transformation rule for connection 



the curvature of trivial connection is obviously trivial: tt(xi, . . . , x^) = 1; 

• the gauge principle does not tell us anything about the evolution of the connection 
itself, so gauge invariant relation describing dynamics of connection {gauge field) 
should be added. 

Let us give two illustrations of how these concepts work in continuous case. 
Electrodynamics. Abelian prototype of all gauge theories. 

Here the set M. is 4-dimensional Minkowski space with points x = (x^) and the set of 
states is Hilbert space of complex scalar (Schrtidinger equation) or spinor (Dirac equation) 
fields tp(x). The symmetry group of the Lagrangians and physical observables is the unitary 
group T = U(l). The elements of T x can be represented as e ~ ta( - x \ 

Let us consider the parallel transport for two closely situated space-time points: 




(43) 



• the curvature of connection 7r(x, y) is defined as the conjugacy clas^3of the holo 
nomy along a cycle of a graph: 



7r(xi, x 2 , . . . , x fc ) = 7r(xi, x 2 )vr(x2, x 3 ) • • • 7r(x fc , xi); 




substituting approximations 



7r(x, x + Ax) 



e 



ip(x,x+Ax) 



1 - iA(x)Ax 



7r'(x, x + Ax) 



e -ip(x,x+Ax) _ 1 _ i A'(x)Ax, 



14' 



The conjugacy equivalence means that ir'(xi, . . . , 



x k ) ~ 7 1 ir(xi, . . . , 



Xk)j for any 7 G V. 
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e - la (x+Ax) _ e -ia( x ) (1 _ iVQ ( x ) Aa; ) ; 
and taking into account commutativity of F = U(l) we obtain 

OOt{ X ) 

A'(x) = A(x) + Va(x) or, in components, A'ix) = AJx) + „ . (44) 

The 1-form yl taking values in the Lie algebra of U(l) and its differential F = (F^) = dA 
are identified with the electromagnetic vector potential and the field strength, respectively. 
To provide the gauge invariance of the equations for field if)(x) we should replace partial by 
covariant derivatives 

<9 M -> Dp = df, - iA^x) 

in those equations. 

Finally, evolution equations for the gauge field A{x) should be added. In the case of 
electromagnetics these are Maxwell's equations: 

dF = first pair (45) 
d ★ F = second pair. (46) 

Here * is the Hodge conjugation {Hodge star operator). Note that equation d46l ) corre- 
sponds to vacuum Maxwell's equations. In the presence of the current J the second pair 
takes the form * d * F = J. Note also that the first pair is essentially a priori statement, 
it reflects simply the fact that F, by definition, is the differential of an exterior form. 

Non-Abelian Gauge Theories in Continuous Space-time. 

Only minor modifications are needed for the case of non-Abelian Lie group F. Again 
expansion of the T-valued parallel transport for two close space-time points x and x + Ax 
with taking into account that tt(x, x) = 1 leads to introduction of a Lie algebra valued 
1-form A = (An) : 

tt(x, x + Ax) w 1 + J 4 At (x)Ax / ". 
Infinitesimal manipulations with formula d43l 

7(x) _1 vr(x,x + Ax)j(x + Ax) — ► t(^)" 1 (1 + A^{x)Ax^) (j(x) + ^gr^ 
lead to the following transformation rule 

A'^x) = 7 (x)- 1 ^(x) 7 (x) + T^)- 1 ^- (47) 

The curvature 2-form 

F = dA + [A A A] 

is interpreted as the physical strength field. In particular, the trivial connection 

*.<*>- 
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is flat, i.e., its curvature F = 0. 

There are different approaches to construct dynamical equations for gauge fields ll25l . 
The most important example is Yang-Mills theory based on the Lagrangian 

L YM = Tr [F A *F] . 

The Yang-Mills equations of motion read 

dF + [A A F] = 0, (48) 
d * F + [A A *F] = 0. (49) 

Here again equation (148T ) is a priori statement — the Bianci identity. Note that Maxwell's 
equations are a special case of Yang-Mills equations. 

It is instructive to see what the Yang-Mills Lagrangian looks like in the discrete approx- 
imation. Replacing the Minkowski space M. by a hypercubic lattice one can see that the 
discrete version of Lym is proportional to ^ . a (7/), where the summation is performed 
over all faces of a hypercubic constituent of the lattice; 

a = 2dimp(r)- x (p(r))-x(y (r)) ; 

where p (T) and p' (F) are fundamental representation of V and its dual, respectively; \ is 
the character; 7^ is the gauge group holonomy around the face /. 

The Yang-Mills theory uses Hodge operation converting fc-forms to (n — A;)-forms in n- 
dimensional space with metric g^ v . In topological applications so-called BF theory plays an 
important role since it does not require a metric. In this theory, an additional dynamical field 
B is introduced. The Lie algebra valued (n — 2)-form B and the 2-form F are combined 
into the Lagrangian Lbf = Tr [B A F] . 



6.2. Quantum Behavior and Gauge Connection 

The Aharonov-Bohm effect (Fig. [9]) is one of the most remarkable illustrations of interplay 
between quantum behavior and gauge connection. Charged particles moving through the 
region containing perfectly shielded thin solenoid produce different interference patterns on 
a screen depending on whether the solenoid is turned on or off. There is no electromagnetic 
force acting on the particles, but working solenoid produces U(l)-connection adding or 
subtracting phases of the particles and thus changing the interference pattern. 

In the discrete time Feynman's path amplitude ll26l is decomposed into the product of 
elements of the fundamental representation p (V) = U(l) of the circle, i.e., of the Lie group 
T = S 1 = R/Z: 

A vm . „p llS) _ ex P (, / U t ) ^ e^o, . . . e ^-M . . . e » W. m 

By the notation Lt-i t we emphasize that the Lagrangian is in fact a function defined on 
pairs of points (graph edges) — this is compatible with physics where the typical La- 
grangians are depend on the first order derivatives. Thus we can interpret the expression 
ir(t - 1, t) = e iLt ~^ £ p (T) = U(l) as U(l)-parallel transport. 
A natural generalization of this is to suppose that: 
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Figure 9. Aharonov-Bohm effect. Magnetic flux is confined within the perfectly shielded 
solenoid; interference pattern is shifted in spite of absence of electromagnetic forces acting 
on the particles. 

• group T may differ from S 1 , 

• dimension of unitary representation p (F) may differ from 1. 

So let us replace expression (l50l ) for Feynman's path amplitude by the following parallel 
transport along the path 



Here a t .t-i are elements of some group T — we shall call it quantizing group — and p is 
an unitary representation of T. Note that in (l50l ) the order of factors is not important due to 
commutativity of U(l). But in (I5TI ) we must use the revers J^l order for consistency with the 
temporal ordering of non-commutative operators. For discrete and especially finite systems 
it is natural to take a finite group as the quantizing group, in this case all manipulations 
— in contrast to the standard quantization — remain within the framework of constructive 
discrete mathematics requiring no more than the ring of algebraic integers (and sometimes 
the quotient field of this ring). On the other hand, the standard quantization can be approx- 
imated by taking 1-dimensional representations of large enough finite groups. 

6.2.1. Illustrative Example Inspired by Free Particle. 

In quantum mechanics — as is clear from the never vanishing expression exp (^S*) for 
the path amplitude — transitions from one to any other state are possible in principle. But 
we shall consider computationally more tractable models with restricted sets of possible 
transitions. 

Let us consider quantization of a free particle moving in one dimension. Such a particle 

■ 2 

is described by the Lagrangian L = m £ . Assuming that there are only transitions to the 
15 This awkwardness stems from the tradition to write operator actions on the left (cf. footnote|4]on page[4]l. 



Ap(r) = P ("t,t-i) • • • P (at,t-i) • • • P (ai,o) • 



(51) 
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closest points in the discretized space we come to the following rule for the one-time-step 
transition amplitudes 

j m{(x+l)~x} 2 -m 
eh 2 = e 4 2/1 

i m(x—x) 2 

eh 2 = l 

j m{(x-l)-x} 2 -m 

eh 2 =e2fi. 

That is, we have evolution rule as an U(l)-valued function R defined on pairs of points 
(graph edges). Symbolically: 

R{x-fx) = 1 G U(l), 
R(x->x-l) = R(x x + 1) = w = e l 2h G U(l). (52) 

Now let us assume that w in (l52l is an element of some representation of a finite group: w = 
p (a) , a £ r = {71 = 1, . . . , 7m}- Rearranging multinomial coefficients — trinomial in 
this concrete case — it is not difficult to write the sum amplitude over all paths from the 
space-time point (0, 0) to the point (x, t) 

t! t ! 

Al(w) = 'S^- / ' x x x — - — -w T . (53) 

Note that x must lie in the limits determined by t: x £ [—£,£]. 

One of the most expressive peculiarities of quantum-mechanical behavior is the destruc- 
tive interference — cancellation of non-zero amplitudes attached to different paths converg- 
ing to the same point. By construction, the sum of amplitudes in our model is a function 
A{w) depending on distribution of sources of the particles, their initial phases, gauge fields 
acting along the paths, restrictions — like, e.g., "slits" — imposed on possible paths, etc. In 
the case of 1-dimensional representation the function A(w) is a polynomial with algebraic 
integer coefficients and w is a root of unity. Thus the condition for destructive interference 
can be expressed by the system of polynomial equations: A(w) = and w M = 1. For 
concreteness let us consider the cyclic group T = Zm = {71, • • • , 7fc, ■ ■ ■ , 7m}- Any of 
its M irreducible representations takes the form p (7^) = w k ~ l , where w is one of the Mth 
roots of unity. For simplicity let w be the primitive root: w = e 2m / AI . 

Fig. [10] shows all possible transitions (with their amplitudes) from the point x in three 
time steps. We see that the polynomial A^ = 3w + 3w 3 = 3w (w 2 + l) contains the 
cyclotomic polynomial Q±(w) = w 2 + l as a factor. The smallest group associated to &i(w) 

— and hence providing the destructive interference — is Z4, which we shall consider as 
quantizing group for the model. 

Fig. [UJshows interference patterns — normalized squared amplitudes ("probabilities") 

— from two sources placed in the positions x = — 4 and x = 4 for 20 time steps. The upper 
and lower graph show interference pattern when sources are in the same (Ac/) = 0) and in 
the opposite (A<fi = tt) phases, respectively. 




x + 1 



>• x 



x — 1 
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Figure 10. Amplitudes for all possible paths in three time steps. 
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Figure 11. Group Z4. Interference from two sources at points -4 and 4. Number of time 
steps T = 20. Phase differences Acp = ^4 — </>_4 between sources are and ir. 

6.2.2. Local Quantum Models on Regular Graphs. 

The above model — with quantum transitions allowed only within the neighborhood of a 
vertex of a 1 -dimensional lattice — can easily be generalized to arbitrary regular graph. 
Our definition of local quantum model on k-valent graph uncludes the following: 

1. Space X = {xi, ■ ■ ■ , x^} is a k— valent graph. 

2. Set of local transitions Ei = {eo,«, e\,i, • • • , e& j} is the set of k adjacent to the vertex 
Xi edges e m>i = (ajj ->■ x m>i ) completed by the edge e 0ji = (xi->Xi). 

3. We assume that the space symmetry group G = Aut (X) acts transitively on the set 

{Ei, • • • , -Ejv x }- 

4. Gioc = Stabc (xi) < G is the stabilizer of Xj. 

5. Oj = {wo,ij " " • j is the se? of orbits of Gi oc on £"j. 

6. Quantizing group T is a finite group: T = {71, • • • , 7a/ }■ 
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7. Evolution rule R is a function on Ei with values in some representation p (T). The 
rule R prescribes p (r)-weights to the one-time-step transitions from Xi to elements 
of the neighborhood of Xj. From the symmetry considerations R must be a function 
on orbits from Oj, i.e., R (e m ^g) = (e mj j) for g £ Gi oc . 

To illustrate these constructions, let us consider the local quantum model on the graph 
of buckyball (see detailed consideration of this graph at page [8]). Here the space X = 
{xi, • • • ,X6o} has the symmetry group G = Aut (X) = Z2 x Alt(5). The set of local 
transitions takes the form Ei = {eo,i, eii, e2.i1 e 3,i}> where 

£•1, i — (^i ^ 3-1, i) ? 

e2,i = (^i -> X2,i) , 
e3,i = (^i -> £3,1) . 

The stabilizer of X{ is Gi oc = Stabc (x;) = Z 2 . The set of orbits of Gi oc on Ei contains 3 
orbits: 

= {uo,i = {eo,i} , = {ei^, e 2 ,j} , w 2j i = {e 3j j}} , 

i.e., the stabilizer does not move the edges (xi — > xi) and (x^ — > xs,i) and 
swaps (xi — > xi t i) and (xj — >■ X2,i) ■ 
The evolution rule takes the form: 

R(xi->Xi) = p(a ), 
R(xi^- x lti ) = R(xi ->■ x 2 ,i) = /o(«i), 
i? (xj -)■ x 3j j) = p(a 2 ), 

where ao 5 «i) "2 £ T. If we take a 1-dimensional representation and move ao — using 
gauge invariance — to the identity element of T, we see that the rule R depends on two 
elements v = p(cti) and w = p{a-2)- Thus the amplitudes in the quantum model on the 
buckyball take the form A(v , w) depending on two roots of unity. To search nontrivial 
quantizing groups one should check (by, e.g., Grobner basis computation) compatibility of 
the system of polynomial equations A(v,w) = &i(v) = &j(w) = 0, where and 
$ j (w) are cyclotomic polynomials. 

6.3. General Discussion of Quantization in Finite Systems 

As is well known, Feynman's approach is equivalent to the traditional matrix formulation 
of quantum mechanics, where the time evolution — > \iPt) of a system from the initial 
state vector to the final is described by the evolution matrix U: = U\ipo). The 

evolution matrix can be represented as the product of matrices corresponding to the single 
time steps: U = Ut^t-i • • • U^-t-i ■ • ■ Ui<-q. In fact, Feynman's quantization — i.e., 
the rules "multiply subsequent events" and "sum up alternative histories" — is simply a 
rephrasing of matrix multiplication. This is clear from the below illustration presenting 
two-time-step evolution of a two-state system (single qubit) in both Feynman's and matrix 
forms — the general case of many time steps and many states can easily be obtained (by 
induction for example). 
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We see that in accordance with Feynman's rules the transition from, e.g., </> 2 to tpi is de- 
termined by the expression 611012 + &i2«22- But this is just the element ui 2 of the matrix 



U\4>), where 



and 



1p2 



product U = B A performing evolution 

Of course, such reduction of sums over histories to matrices is applicable to the case of 
transitions along paths being gauge connections as in d5TT ). In this case matrix elements of 
an Nx x Nx evolution matrix U are themselves matrices from the representation p (T). We 
can ignore this particular block structure of the matrix and consider U as an N x N matrix 
over the field C, where N = Nx x dim p (F). 

In quantum mechanics, the evolution matrices U are unitary operators acting in Hilbert 
spaces of state vectors (called also "wave functions", "amplitudes" etc.)- Quantum mechan- 
ical particles are associated with unitary representations of some groups. According to their 
dimensions, these representations are called "singlets", "doublets", etc. Multidimensional 
representations describe the spin. A quantum mechanical experiment is reduced to com- 
parison of the system state vector | ip) with some sample state vector | <p) . According to the 
Born rule, the probability to observe coincidence of the states is equal to |(0|^)| 2 , where 
(•|-) is the inner product in the Hilbert space. To see what these constructions may look like 
in the constructive finite background, let us assume that evolution operators are elements of 
a representation of a finite group. 



6.3.1. Permutations and Linear Representations 

Having a finite group G = {e\ = 1, . . . , e m }, we can easily describe all its transitive ac- 
tions on finite sets [27]. Any such set = {u}\, . . . ,to n } is in one-to-one correspondence 
with the right (or left) cosets of some subgroup H < G, i.e., U = H \ G (or G/H) is the 
homogeneous space (or G-space). Action of G on Q, is faithful if the subgroup H does not 
contain normal sugroups of G. We can write actions in terms of permutations 

Maximum transitive set Q, is the group itself, i.e., in the above construction H = {1}. The 
action on Q = G is called regular and can be represented by permutations 



n (s) = ( &% )> i = l,...,m. 



(54) 
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To introduce "numerical" ("statistical") description, let us assume that oVs are basis ele- 
ments of a linear vector space % over a field F 



H = Span (wi, ■ 



(55) 



i.e., we prescribe F- valued "weights" to the elements oji € O. Then we can write permuta- 
tions in the matrix form: 



7T 



(s) -> P(g) = [p(g)ij] , where p{g\ 



•J 



1, 



,n; 



(56) 



1, if q = /3, 



for a,j3 £ tl. 



[0, if (3 

The so defined function p is called a permutation representation. The matrix form of ((54 



n( 5 ) -> P( 5 ) = [P(^] , Pfojy = <y, 



1, 



, m 



(57) 



is called the regular representation. It is assumed that F is an algebraically closed field — 
usually the field of complex numbers C. But in the case of finite groups the quotient field 
of the ring A of algebraic integer!^ fT2| is sufficient for all reasonable purposes — A is a 
constructive subset of C. 

Let us recall some relevant background information about linear representations of finite 
groups Il28l . 

1 . Any linear representation of a finite group G is unitary since there is always an unique 
invariant inner product (-|-) making any space of representation % into a Hilbert 
space. 

2. All possible irreducible unitary representations of the group G are contained in the 
regular representation (157T ). More specifically, all matrices ((571 ) can simultaneously 
be reduced by some unitary transformation S to the form 



A1G7) 



fA 2 (fl) 



A2G7) 



d r < 



(58) 



Here r is the number of different irreducible representations Aj of the group G. This 
number coincides with the number of conjugacy classes^ in G. The number dj is 



16 The ring of algebraic integers consists of the roots of monic polynomials with integer coefficients. A 
polynomial is called monic if its leading coefficient is unit. 

17 The jth conjugacy class Cj C G consists of all group elements of the form g~ 1 Cjg, where Cj £ Cj is 



some (arbitrary) representative of the class, g £ G, j 



,,r. 
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simultaneously the dimension of Aj and its multiplicity in the regular representation, 

so it is obvious that df + d^-l + d\ = \G\ = m. It can be proved also that any dj 

divides the number of elements of G: dj \ m. 

3. Any irreducible representation Aj is determined uniquely (up to isomorphism) by its 
character Xj ■ The character is a function on G defined as Xj ( a ) = Tr A j (a), a G 
G. The character is a central or class function, i.e., it is constant on the conjugacy 
classes: Xj( a ) = Xj(9~ la 9): a,g £ G. Any class function tp(a) on G is a linear 
combination of the characters xii ■ ■ ■ > Xr- 

4. All values of Xj and eigenvalues of Aj are elements of the ring A of algebraic inte- 
gers, moreover the eigenvalues are roots of unity. 

5. A convenient form of describing all irreducible representation of a finite group G is 
the character table. The columns of this table correspond to the conjugacy classes 
of G while its rows correspond to the characters Xj of the inequivalent irreducible 
representations of G. 





1 c 2 


C r 


Xi 


1 1 


1 


X2 


X2(ci) X2(C 2 ) • 


•• X2(c r ) 


Xr 


Xr(Cl) Xr{c 2 ) ■ 


■■ Xr(c r ) 



The jth column is indicated by a representative Cj G Cj of the jth conjugacy class 
Cj. Conventionally we take c\ = 1 and xi to be the trivial character corresponding 
to the trivial 1 -dimensional representation. 

6.3.2. Interpretation of Quantum Description in Finite Background 

Let us discuss sketchy (more detailed presentation see in |29l ) constructive approach to the 
interpretation of quantum description. 

Summarizing the above, we see that dynamics of finite quantum model of any type is 
reduced ultimately to a single finite-dimensional unitary k x k matrix U describing transi- 
tions between initial and final vectors in some /c-dimensional Hilbert space In the finite 
background the matrix U is an element of unitary representation A of a finite group G, i.e., 
the number of all possible evolutions is equal torn = \G\. We shall assume, as is accepted 
in quantum mechanics, that A is direct sum of irreducible representations Aj from (I58T ). 
The decomposition of the Hilbert space into irreducible components is an important part of 
the mathematical formulation of quantum mechanics. Such dependence on the choice of 
the basis in the Hilbert space may seem unusual for a physical theory. But, in fact, a basis 
in which the Hilbet space is reduced — we shall call such a basis quantum basis — simply 
reflects the structure of underlying symmetry group. 

We can construct a G-space Vt = {ui, . . . ,o; n ]0 in such a way that its permutation 
representation (1561 ) contains A as subrepresentation (obviously n > k). That is, the space 

18 In the case that A is reducible representation, the set Q, may be intransitive union of transitive G-spaces. 
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Tik is subspace of the Hilbert space % n of the permutation representation. We shall call the 
basis {loi, . . . , u n } in the space % n the permutation basis. Transitions from the permutation 

to quantum basis for matrices U and vectors ip)£ Un are given by the formulas 



S^UpS, 
1 



s 



(59) 
(60) 



Now we can embed any evolution U with the matrix A in the space Hk into the evolution 
U in the space % n . In the quantum basis the matrix of U takes the form 



U„ 



A 0' 
A 



(61) 



where A is an {n — k) x (n — k) matrix. Due to the form of (|6"TT) the evolution U de- 
scribed by A is completely independent of the components of vectors of % n related to A. 
The "hidden variables" that can come from the additional components describe degrees of 
freedom reflecting indistinguishability of u>i's lying in the same group orbit. The evolution 
U is simply a permutation of cjj's and can not manifest anything quantum. 



Illustration. A quantum model with the group Sym(3). The group G = Sym(3) 
is the group of all permutations of three objects. This is the smallest non-commutative 
group. Its 6 elements form the following 3 conjugate classes C\ = {1 = ()}, C2 = 
{01 = (12), a 2 = (23), 03 = (13)}, C 3 = {h = (123), b 2 = (132)}. We used here the 
cyclic notation for permutations. The group has three nonequivalent irreducible representa- 
tions described by the character table 





1 


a,i 


h 


Xi 


1 


1 


1 


X2 


1 


-1 


1 


X3 


2 





-1 



Let us take for example the 2-dimensional representation A with the character x,3- The 
representation is given by the following set of 2 x 2 matrices: 

"l 0" 
1 



A(l) 



A(ai) 



e~ 27ri / 3 



A(6i 



3 27ri/3 




, A(a 2 


-2ni/3 



1 

1 



, A(6 2 ) 



A(o 3 ) = 

e -27ri/3 





e 2 "/ 3 

e -2ni/3 Q 



e 27ri/3 



The regular permutation representation of Sym (3) is 6-dimensional. But 3-dimensional 
faithful permutation representation induced by the action on the homogeneous space 
Sym (2) \ Sym (3) = Q, = {011,012,003} also contains A. Since any permutation repre- 
sentation contains trivial 1 -dimension subrepresentation, the only possible choice of the 
addition A is the representation corresponding to the first row of the above character table. 
Thus, for U q we have 
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U q (l) 



1 
1 
1 



UJ ai ) 



e" 2 ™/ 3 0' 
e 27ri/3 o o 

1 



, U q (a 2 ) 



1 

1 
1 



i U q (a 3 ) 



e 2 ™/ 3 0' 



e -27ri/3 





u q (h) 



e 27Ti/3 








e -27Ti/3 Q 

1 



u q (b 2 ) 



■ e -2iri/3 q 0' 

e 2 ™/ 3 
1 



In the permutation basis we have 



U P (1) 



1 
1 
1 



U p (ai) 



1 

1 
1 



U p (a 2 ) 



1 
1 
1 



U p (a 3 ) 



1' 

10 

1 



1 

1 

1 



U p (b 2j 



1 

1 
1 



The most general unitary matrix of transition from the permutation to the quantum basis 
takes the form 



i 

e 2vri/3 



1 



-27ri/3 Ri 



-2ni/3 „27ri/3 



J8<- 



, where a, /3 are arbitrary real parameters. (62) 



Any quantum evolution of the form 



U \ 4>), where 



and U is one of the matrices A; can be extended to the evolution 



and 



U a 



_1p2 

■j q ), where 



and 



■03 



, 03 is arbitrary additional component. Then, applying the 



transformation S, we come to the classical evolution with the matrix U p = SU q S 1 which 
simply permutes the components of the initial vector 



S 



7! 



01 + 4> 2 + e^3 

e 2^/3^ 1 + e -2vri/3^ 2 + ^ife 
e -2 7 ri/30 l + e 27ri/3^ 2 + Qpife 



without performing any algebraic manipulations with the components. 
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7. Conclusion 

In this chapter we discuss the general concept of discrete dynamical system and its special- 
ization involving underlying space structures. We apply various constructive approaches to 
study discrete and finite dynamical systems. 

We construct a family of groups unifying space and internal symmetries in a natural 
way. This construction generalizes the standard direct and wreath products. 

We introduce the concept of a system of discrete relations on an abstract simplicial 
complex. This system can be treated as a natural generalization of cellular automata or as a 
set-theoretical analog of systems of polynomial equations. 

We developed and implemented algorithms for analyzing compatibility of systems of 
discrete relations and for constructing canonical decompositions of discrete relations. 

Applying the technique described above to some cellular automata — a particular case 
of discrete relations — we obtained a number of results. The most interesting among them, 
in our opinion, is the demonstration of how the presence of non-trivial proper consequences 
may determine the global behavior of an automaton. 

We suggest an algorithmic approach — based on discrete symmetry analysis and im- 
plemented in C — for construction and investigation of discrete dynamical models — de- 
terministic, mesoscopic and quantum. We hope that our approach can be used in various 
practical applications, such as, for example, simulation of nanostructures with nontrivial 
symmetry properties. 

We demonstrate that soliton-like moving structures — like "spaceships" in cellular au- 
tomata, "traveling waves" in mathematical physics and "generalized coherent states" in 
quantum physics — arise inevitably in deterministic dynamical systems whose symmetry 
group splits the set of states into finite number of group orbits. 

We formulate the gauge principle in the form most suitable for discrete and finite sys- 
tems. We also propose a method — based on introduction of unitary gauge connection of 
a special kind — for quantizing discrete systems and construct simple models for studying 
properties of suggested quantization. 

We show that if unitary operators describing dynamics of finite quantum system form 
finite group, then the system can be embedded into a classical system with a simple behav- 
ior. We hope that discrete and finite background allowing comprehensive study may lead 
to deeper understanding of the quantum behavior and its connection with symmetries of 
systems. 

To study more complicated models we are developing C programs based on computer 
algebra and computational group theory methods. 
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